Concurrence and entanglement entropy in a dimerised spin-1/2 two-leg ladder 
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We consider isotropic spin-1/2 two- leg ladders with dominant spatially-modulated rung exchanges. 
We study the effect of a uniform magnetic field on the ground state phase diagram of the model 
using perturbation theory and the numerical Lanczos method. The ground state phase diagram 
consists of two gapless Luttinger liquid (LL) and three gapped phases. Numerically, we calculate the 
concurrence between two spins and the entanglement entropy between legs. Numerical experiment 
shows that the gapless LL phases are fundamentally different. In the first LL phase, only spins 
on rungs are entangled, but in the second LL phase the spins on legs are long-distance entangled. 
Therefore, the concurrence between spins on legs can be considered as a function to distinguish the 
LL phases. 

PACS numbers: 75.10.Jm; 75.10.Pq 



I. INTRODUCTION 

Entanglement is a quantum physical property that de- 
scribes a correlation between quantum systems that has 
no classical analogy. When a quantum system is entan- 
gled, the best possible knowledge of the whole does not 
include the best possible knowledge of its parts. In par- 
ticular, a quantum system possesses additional correla- 
tions that do not have a classical counterpart. The study 
of the nature of these quantum correlations has attracted 
a lot of interest in the last decade. 

At zero temperature, quantum fluctuations play the 
dominant role and with the ground state properties the 
physics of a system determined . The quantum fluctua- 
tions cause a fundamental change in the state of a system 
which is known as the quantum phase transition. Entan- 
glement is expected to play an essential role in quantum 
phase transitions, where quantum fluctuations manifest 
themselves at all length scales. It is believed that the 
entanglement estimators give important insight into the 
physics of the quantum phase transitions. 

Low-dimensional quantum spin systems exhibit quan- 
tum phase transitions for different driving parameters, 
like a magnetic field. In particular, we address this 
problem by selecting a system of spins which are ar- 
ranged on a two-leg ladder lattice with antiferromag- 
netic isotropic couplings along the legs and rungs. The 

= ^ isotropic antiferromagnetic two-leg ladders, have 
attracted much interest for a number of reasons. On 
the one hand, there has been remarkable progress in re- 
cent years in the fabrication of such ladder compounds^. 
On the other hand, spin-ladder systems pose interesting 
theoretical problems, since they have a gap in their ex- 
citation spectrum. In the presence of a magnetic field, 
they reveal an extremely rich behaviour, dominated by 
quantum effects. These quantum phase transitions have 
been investigated intensively both theoretically^r— and 
exp er iment all jsi^"— . 

In a recent paper, the ground state phase diagram of a 
different kind of two-leg ladder with spatially-modulated 
spin-exchange coupling constants has been studied^^ in 



the limit of strong rung exchange and magnetic field. The 
Hamiltonian of the model is defined as 

n , a n . Q 

+ J'±5I[l + (-l)"'^]S"d •Sn,2, (1) 



where Sn.Q is the spin S 
(n=l,...,L) and leg a [a — 
is antiferromagnetic, 



= 1/2 operator on rung n 
1,2). The interleg coupling 
Jj_(l ± 5) > 0. The non- 



modulated case^^i^i, (5 = 0, undergos two quantum phase 
transitions at critical fields, ffci and . The first cor- 
responds to the transition from the gapped rung-singlet 
phase to the gapless LL phase. The second represents the 
transition from the LL phase into the fully polarised fer- 
romagnetic phase. Recently, the bipartite entanglement 
measures close to the quantum critical points^^ have been 
investigated. It was found that the derivatives of the en- 
tanglement measures diverge as H — -ffc2 but remain 
finite as H — > H^^ . The entanglement entropy was also 
studied for two-leg ladders within spin wave theory and 
the DMRG method^^. 

For the modulated case, (5 7^ 0, the continuum-limit 
bosonizatin analysis^^ shows two additional quantum 
phase transitions between a middle gapped phase and 
two gapless LL phases. These transitions manifest them- 
selves most clearly in the presence of a new magnetisation 
plateau at magnetisation equal to one half of its satu- 
ration value. Using the numerical Lanczos method, the 
existence of four critical fields, the magnetisation plateau 
at half of the saturation value and two gapless LL phases 
are confirmed^i. 

In this paper we continue the study of this model. 
First, using perturbation theory, we find four critical 
fields in good agreement with the previous results. Then, 
using the numerical Lanczos method, the concurrence 
and entanglement entropy will be calculated and used 
to show that the LL phases are fundamentally different. 

The outline of the paper is as follows. In the forth- 
coming section we present the perturbation results. In 
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section III we present the results of an accurate numeri- 
cal experiment on the entanglement between spins at zero 
temperature. Finally, we discuss and summarise our re- 
sults in section IV. 



II. PERTURBATION RESULTS 

We are interested in the behaviour of the model Eq.([T]) 
in the limit of strong rung exchange, J±,H 3> J\\. For 
this aim, it is convenient to rewrite the Hamiltonian given 
in Eq.([T]) in the form 



the gapped rung singlet phase to the LL phase. The 
second one. He., represents the transition from the LL 
phase to the fully polarised ferromagnetic phase^^iSi. 

In the case of dimerised ladder o^^i^"^ , 5^0, there are 
two kind of rungs: strong and weak. In this case, by 
increasing the magnetic field, first weak rungs start to 
melt and undergo into the triplet state with respect to 
the strong rungs. Therefore, it is natural to find two ad- 
ditional critical fields in compare the non-dimerised case 
{S = 0). However, we found that second order pertur- 
bation theory should be used for solving the eigenvalue 
equation. The obtained critical fields are 



V, 



-Ho = Ji^[l + (-l)"<S]S„,i-S„,2 

n 
71, a. 

^ ^ "^ll ^ ^ ^n,a ■ ^n+l.ct- 



(2) 



The unperturbed part, "Ho, is the Hamiltonian of a non- 
interacting N/2 pair of spins. The eigenstate of the un- 
perturbed Hamiltonian is written as a product of pair 
states. By solving the eigenvalue equation of an indi- 
vidual pair, one can easily find the eigenstates as \S) = 
^(1 it) - I n)), iTi) = I tt), |To) = ^(1 it) + I ti)), 
|T_i) = I it). Their eigenvalues are, respectively. 
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An) 



\Ti) - jJ±(''^)~-" 7 ^\To} ~ 4' 

H. When the magnetic field is zero. 



the ground state of a pair of spins is singlet. By increas- 
ing the magnetic field, the energy of jTi) becomes closer 
to E^s) and an abrupt quantum phase transitioni occurs 
at He = J±{n). By further increasing the field, |Ti) will 
become the ground state of a distinct pair. Thus, for 
a strong enough magnetic field, we have a situation in 
which the behaviour of the system is determined by the 
singlet IS") and triplet |Ti) states. 

Let us start with the case of (5 = 0. Since the ground 
state energy of a distinct pair is twofold degenerate at 
He = Jl, the ground state energy of the unperturbed 
Hamiltonian, Hq, is 2^ times degenerate^. The pertur- 
bation V splits this degeneracy. By applying first or- 
der perturbation theory for finite ladders (L = N/2 ~ 
2,3,4,5,6) with periodic boundary conditions and gen- 
eralising the results to the thermodynamic limit we have 
found the critical fields to be 



He, = J±+2J||, 



(3) 



which arc in good agreement with the results obtained 
using analytical field theor y^^i^^ . It was shown that the 
first critical field, H^ , corresponds to the transition from 



He, 


= (1 


Hc- 


= (1 


He+ 


= (1 


Hc2 


- (1 



= (1-<5)J±- 



2SJj_' 



(l + <5)Jj_ + Jii + 



2(5 J I 



(4) 



We emphasise that the above perturbation results 
are in good agreement with estimations made within 
the continuum-limit approach^^ and numerical Lanczos 
method^!. It has been show n^^i^^ that the ground state 
phase diagram consists of five phases: (I.) the rung- 
singlet phase in the region H < H^, (II.) the first LL 
phase in the region H^ < H < H^- , (HI.) the gapped 
mid-plateau state in the region H^- < H < i/c+- (IV.) 
the second LL phase in the region He+ < H < Hc2 
and (V.) the saturated ferromagnetic phase in the region 
H > He,. 

On the other hand, we also calculated the ground state 
of the system in different sectors of the ground state 
phase diagram, which allows us to follow the melting 
process on the entanglement function^i^. The entan- 
glement of formation is defined as 



E = -Xlog2X - (1 - X)log2{l - X), 



(5) 



where X = ^(1 + vT— C^) and C is the concurrence 
which is given by 



Ci„. = 2max{0,ci^lcl'„l}, 



where 



^(1) 




-9Zf 












r(2) _ 




+9zr 


+ (.9;,n - 9Z) 



L)2 



(6) 



9i!n ~ {^?^m) denotes the correlation function between 
spins on sites / and m and = {S[). The notation (...) 
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represents the ground state expectation value. Because 
J± ^ J|| , the entanglement of formation between spins 
on a rung is an important quantity. In the region of the 
magnetic filed, H < H,.^, the ground state is written as 
a product of singlet rung states. In this region the con- 
currence between spins on a rung attains its maximum 
value, one. For very large values of the field, H > Hc^, 
the system is in the saturated ferromagnetic phase and 
the ground state is a product of triplet (|Ti)) rung states. 
In the saturated phase the concurrence takes zero value. 
Therefore, during the melting process, the concurrence 
between spins on rungs changes from the completely en- 
tangled state into the non-entangled state. In addition, 
using the perturbation results on finite ladders, we found 
a relation for the concurrence between spins on a rung in 
different subspaces versus the number of rungs as 

where i is the number of triplet rungs in a subspace. In 
the rung singlet phase, the number of triplet rungs is zero 
and the concurrence is maximum. On the other hand, in 
the saturated ferromagnetic phase all rungs are triplet, 
i = L, and particles are not entangled (C = 0). 
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FIG. 1: (colour online). The concurrence between spins on 
rungs as a function of the external magnetic field for ladders 
with lengths TV = 20, 24, 28 and uniform rung couplings, 5 = 
0, J± = 6.0 and J|| = 1.0. 



III. NUMERICAL EXPERIMENT 

By performing an experiment, one can obtain a clear 
picture of the entanglement phenomena in the ground 
state magnetic phases of the model. Since a real experi- 
ment cannot be done at zero temperature, the best way 
is to undertake a virtual numerical experiment. A very 
famous and accurate method in field of numerical ex- 
periments is known as the Lanczos method^i^. In this 
section, to explore the nature of the spectrum and the 
quantum phase transition, we used the Lanczos method 
to numerically diagonalise ladders with lengths up to 
N = 2L = 28 and different values of the exchanges^l. 
The energies of the lowest-energy eigenstates were ob- 
tained for ladders with periodic boundary conditions. 



We have also considered the case of a ladder with alter- 
nating rungs S 0. Critical fields calculated numerically 
for ladders with exchange parameters J± = -y, S = jj, 
J|l = 1 are22 

He, = 4.47 ±0.01, 
He- = 5.34 ±0.01, 
He+ = 6.87 ±0.01, 

= 7.62 ±0.01, (8) 

which are in good agreement with the second-order per- 
turbation results (Eq.Q). 

B. CONCURRENCE 



A. CRITICAL FIELDS 

Since in a critical field the energy gap should be closed, 
the best way to find the critical fields is to investigate the 
energy gap, which is recognised as the difference between 
the energies of the first exited state and the ground state 
in finite two-leg ladders. As previous section, we start 
our consideration with the case of isotropic and uniform 
rung exchange ladders {5 = 0). We have calculated the 
energy gap for ladders with exchanges J\\ = 1, J± = 
6 and lengths N = 12, 16, 20, 24, 28. By means of the 
phcnomcnological rcnormalisation group tcchnique^^ we 
have found H^i = 5.04 ± 0.01 and = 8.12 ± 0.01, in 
good agreement with the perturbation results (Eq.([3])). 



In this section we focus on the entanglement of for- 
mation as a measure of the entanglement. We compute 
the entanglement between two spins, which is known as 
the concurrence (Eq.®). In the region H < H^, where 
the system is in the non-magnetic rung-singlet phase, the 
correlations (7^^, g^^, g^^ should be the same and take the 
value —1/4, and therefore the concurrence must take the 
maximum value C = 1. On the other hand, for very large 
values of the magnetic field, H > Hc^ , the ground state 
of the system is written as 

|Gs) = int-t). (9) 

In this saturated ferromagnetic phase, the correlations 
arc g^^^ = = 0,5^^ = 1/4 and = 1/2. Therefore, 
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FIG. 2; (colour online). The concurrence between spins on 
strong (SR) and weak (WR) rungs versus applied magnetic 
field for dimerised ladders with lengths A'' = 16, 24 and ex- 
change parameters ^ = XT> "^x = a^nd Jy = 1.0. 



FIG. 3: (colour online). The concurrence between nearest 
neighbour (NN) spins on legs as a function of the magnetic 
field H for dimerised ladders with lengths N = 16, 20, 24, 28 
and exchange parameters ^ = XTi -^x = a-nd Jy = 1.0. 



the concurrence takes the minimum value C = 0. 

The numerical Lanczos results describing the concur- 
rence for a two-leg ladder with antiferromagnetic legs and 
the same rungs {5 = 0) are shown in Fig.[TJ In this figure 
the concurrence between two spins on a rung is plotted 
as a function of H for ladder lengths N = 20, 24, 28 and 
exchanges Jy = 1, Jl = 6. It can be seen that the nu- 
merical results are in reasonable agreement with what is 
expected. Essentially our numerical experiment shows 
that in the absence of a magnetic field, spins on a rung 
are maximally entangled. By applying a magnetic filed 
up to the first critical field Hc^ , the concurrence remains 
constant. This behaviour is in agreement with expecta- 
tions, based on the general statement that in the gapped 
rung-singlet phase, the change in any physical function 
appears only at a finite critical value of the magnetic 
field equal to the spin gap. The concurrence starts to de- 
crease as soon as the field passes beyond the first critical 
field iJci . Indeed the quantum correlations of the two 
spins on rungs decrease with increasing magnetic field. 
In the intermediate region the decreasing behaviour of 
the concurrence continues up to the second critical field 
, where it takes the zero value. As expected in the re- 
gion of the saturated ferromagnetic phase, H > , the 
concurrence remains zero. The oscillations of the con- 
currence at finite N in the region 7?^ < H < Hc2, arc 
the result of level crossings between the ground state and 
excited states of the model. 

Let us now present our numerical results for the 
dimerised case 6^0. In Fig. [5] we have plotted the 
concurrence of two spins on strong ( J_l(1 + 5)) and weak 



{J±{1 — S)) rungs as a function of the magnetic field H. 
We have considered two-leg ladders with different lengths 

= 16,24 and exchange parameters J± = S = jj, 
Jy = 1. It can be seen that spins on weak and strong 
rungs are maximally entangled at iJ = 0. By increas- 
ing the magnetic field from zero, the spins on strong and 
weak rungs remain maximally entangled up to the first 
critical field H = ■ For H > the concurrence 
between spins on weak rungs drops very rapidly when 
compared to the concurrence between spins on strong 
rungs. Indeed, the quantum correlations of two spins 
on strong rungs and weak rungs decrease with increas- 
ing magnetic field, but with the different intensity. In 
the intermediate gapless LL region < H < H~, the 
quantum correlations of strong and weak rungs diminish 
down to H = H~ and the concurrences reduce to ~ 0.8 
and ^ 0.15 respectively. In the half-plateau state, the 
gap of the system is rc-opcned and a plateau emerges 
in the curve of concurrences with the values 0.85 and 
~' 0.15. By further increasing the magnetic field and for 
H > , the concurrence between spins on strong rungs 
drops very rapidly when compared to the concurrence on 
weak rungs. Finally, in the full-saturated ferromagnetic 
state, all of concurrences disappear and the entanglement 
of the state is exactly zero. 

Apart from the concurrence between spins on rungs, it 
is completely natural to ask about the concurrence be- 
tween spins on legs. To find a clear answer, we have cal- 
culated the concurrence between spins on the legs. We 
have plotted the concurrence between two spins which 
are nearest neighbours (NN) on a leg in Fig. [3] The nu- 
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FIG. 4: (colour online). The concurrence between next 
nearest neighbour (NNN) spins on legs as a function of the 
magnetic field H for dimerised ladders with lengths A*' = 
16,20,24,28 and exchange parameters 5 = jj, Jx = ^ and 
J,i = 1.0. 



FIG. 5; (colour online). The concurrence between spins on 
legs in the second LL phase versus the separation distance 
n for a ladder with length N — 28 and exchanges 5 = j^-, 
Jx = ^ and J|| = 1.0. 



merical results in this figure are computed for two-leg 
ladders with different lengths TV = 16, 20, 24, 28 and ex- 
change parameters J± = 6 = jj, J\\ =1. It is clearly 
seen that the NN spins on legs are not entangled in the 
absence of a magnetic filed. By applying a magnetic field, 
the concurrence between NN spins remains zero up to the 
critical field H = . By increasing the magnetic field 
beyond this value, the NN spins on legs become entan- 
gled in the plateau state H^- < H < 11^+. This means 
that the magnetic field increases the quantum correla- 
tions of the two spins which are NN on a leg. This is a 
dual effect of the magnetic field, in which increases in the 
quantum correlations of two NN spins on legs accompany 
decreases in the quantum correlations of spins on rungs. 
Indeed, in the plateau state there are three types of quan- 
tum correlations in the system. These correlators are the 
source of the mid-plateaus in the different parameters of 
the system, such as magnetisation, with further increases 
in the magnetic field H > , the concurrence between 
NN spins decreases and becomes zero at the saturation 
magnetic field H = ■ 

To get more information about the entanglement be- 
tween spins, we have also calculated the concurrence be- 
tween next nearest neighbours (NNN) on legs. The nu- 
merical results on the NNN spins on a leg are plotted in 
Fig. H for two-leg ladders with lengths N = 16, 20, 24, 28 
and exchange parameters Jx = (5 = jr, Jy = 1. It 
is completely clear that the NNN spins are not entan- 
gled at _ff = 0. In the presence of a magnetic field the 
NNN spins will be entangled only in the second gapless 



LL phase, 11^+ < H < H^.^- This is related to the fact 
that in the second LL phase, the magnetic field increases 
the quantum correlations of the two spins which are NNN 
spins on a leg. The concurrence between them is equal 
to zero in the mid-plateau state. Our numerical results 
show that only NN spins on legs are entangled in the 
mid-plateau state. As in the case of NN spins, the con- 
currence between NNN spins decreases and becomes zero 
at the saturation magnetic field H — ■ To gain fur- 
ther insight into the second LL state, we have computed 
numerically the entanglement between the spins with dif- 
ferent separation distances for ladders with lengths up to 
N = 28. The concurrence between an arbitrary spin and 
another spin, say Sn has been plotted versus the sepa- 
ration distance, n, in Fig. [51 It is very interesting that 
the entanglement of two such spins on legs does not van- 
ishes by increasing the distance between them. It means 
that, in the second LL state, < H < Hc2, the range 
of the quantum correlations between spins on the legs is 
long-distance. There is a non- vanishing long-distance en- 
tanglement between two spins on legs in dimerised ladder 
systems. The oscillation with period two is related to the 
fact that in the dimerised ladders there are two kind of 
rungs. 



C. ENTANGLEMENT ENTROPY 

Up to this point we have presented the results of the 
numerical calculations on the entanglement of formation, 
which is useful for studying the entanglement between a 
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pair of particles. However in some cases, two groups of 
particles can be entangled. Measurements of entangle- 
ment between subsystems, chiefly using entropic quanti- 
ties, have an advantage over traditional two-point corre- 
lation functions in that they encode the total amount of 
information shared between the subsystems without the 
possibility of missing hidden correlations. 

The entanglement entropy (EE) is defined as the von 
Neumann entropy of a reduced density matrix of a sub- 
system. This function is useful for understanding the 
different quantum phases in condensed matter systems. 
The EE is defined as 

E-^ ^-{\ogpA) = -TrA[pA\ogpA\, (10) 

where, pA = TrB[p\ and p is the density matrix of the 
ground state \Gs). It is assumed that the system consists 
of subsystems A and B. The EE quantifies the informa- 
tion describing how much the state is entangled between 
the subsystems A and B. 

In our numerical experiment, legs are selected as sub- 
systems A and B. The Lanczos numerical results on the 
function E""^ , are plotted in Fig.lSJa) for two-leg ladders 
with lengths N = 12, 16, 20 and exchange parameters 
Jj_ = (5 = jj, J|| = 1. It is Clearly seen that the 
legs are entangled in the absence of an external magnetic 
field. They remain entangled in the presence of a mag- 
netic field up to the first critical field H^- It has been 
found that, away from the critical points, the EE scales 
with the size of the boundary between A and B, which 
is known as the area law^. It is completely clear that 
in the rung-singlet phase, H < Hci, the EE is propor- 
tional to the number of rungs, which is known as the area 
law in ladder sytems^i^. On the other hand, a different 
concept of entanglement has been suggested, known as 
valence bond entanglement entropy^i^, E^^ = ln(2).p, 
where p is the number of singlets crossing the bound- 
ary between regions A and B. It is also clearly seen 
that in the singlet region, H < , the valence bond 
entropy is less than the von Neumann entropy, in good 
agreement with Ref.f^]. As soon as the magnetic field 
increases from the first critical field, H > Hc^, the en- 
tanglement entropy decreases and the process stops at 
the critical field Hf,^ . Like the other physical functions, 
a plateau appears in curve of the EE function in the in- 
termediate gapped state, H^- < H < Hf.+ . Finally, by 
further increasing the magnetic field, H > Hc+, the EE 
decreases and vanishes at the saturation magnetic field 
Hc2- The area law suggests that one can find a better 
picture by dividing the entanglement entropy, E"^ , by 
the number of rungs, L. In Fig. we have therefore 
plotted E'"^ /L as a function of the magnetic field H. It 
is clear that legs are maximally entangled in the region 
H < , and that in the intermediate gapped plateau 
region, H^- < H < i?c+j the entanglement entropy 
is almost equal to half of the saturation value. There 
is clearly a decreasing process in both the LL regions 
He, <H < He- and Hc+ < H < Ha^. 
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FIG. 6: (colour online.) (a) The entanglement entropy 
between two legs for dimerised ladder with lengths N = 
12,16,20 and exchanges 5 — -j, Jx = ^ and Jy = 1.0. 
(b) The entanglement entropy per rungs as a function of the 
magnetic field H. 



IV. CONCLUSION 

In this paper, we have focused on the entanglement 
between spins in isotropic spin- 1/2 two- leg ladders with 
dominant spatially-modulated rung exchanges. Using 
perturbation theory and the numerical Lanczos method 
we have studied the effects of an external magnetic field 
and a space modulation on the concurrence and the en- 
tanglement entropy. By employing perturbation theory. 



7 



four critical fields in good agreement with the previous 
Bosonization results are found. Using the perturbation 
results on finite ladders, we have found a relationship be- 
tween the concurrence between spins on a rung and the 
number of rungs (Eql?]). 

Moreover, we implemented the Lanczos method to nu- 
merically diagonalisc ladders with finite length up to 
N = 28. Using the exact diagonalisation technique, wc 
calculated the energy gap, concurrence and the entangle- 
ment entropy between legs for different values of the ex- 
ternal magnetic field. In the intermediate gapped state, 
we found that a non-zero plateau also appears in the 
plot of concurrence between spins on rungs and nearest 
neighbour spins on legs. On the other hand, our numer- 
ical results on the concurrence showed that the two gap- 
less LL phases in the ground state phase diagram of the 



model are fundamentally different. In one of them only 
spins on rungs are entangled, but in the next LL phase 
the spins on legs are long-distance entangled. Thus the 
concurrence between spins on legs can be considered as 
a function to distinguish the LL phases. 
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